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Abstract. We analyze multipartite entanglement between atomic ensembles within 
quantum matter-light interfaces. In our proposal, a polarized light beam crosses 
sequentially several polarized atomic ensembles impinging on each of them at a 
given angle on. These angles are crucial parameters for shaping the entanglement 
since they are directly connected to the appropriate combinations of the collective 
atomic spins that are squeezed. We exploit such scheme to go beyond the pure state 
paradigm proposing realistic experimental settings to address multipartite mixed state 
entanglement in continuous variables. 
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1. Introduction 

Experiments exploiting the quantum character of atom-light interactions are often 
carried out in the regime of cavity QED where, due to cavity, the coupling between 
individual atoms and photons, represented by the optical depth is highly enhanced. An 
alternative approach to achieve an efficient atom-photon coupling are atomic ensembles, 
where large optical depths are obtained due to to the macroscopic or mesoscopic number 
of atoms of the sample, as proposed initially by Kuzmich and coworkers pQ. 

Atomic ensembles typically refer to samples composed of 10 6 to 10 12 atoms. In 
the atomic ensembles the spin of each atom, j i; originating either from its hyperfme 
structure or from its nuclear spin (e.g. alkaline-earth fermions) can be disregarded and 
the sample can be characterized by its collective spin J = (J x , J y , J z ) = £\ jj. If the 
sample is polarized, the component of J along the polarization axis can be approximated 
to a c-number, let's say J x = Y^iJx — (Jx), while the orthogonal components behave as 
conjugate canonical variables [J y , J z ) = ih(J x ). Such an elegant and simple description of 
polarized atomic ensembles makes them extremely appealing for quantum information 
processing in the domain of continuous variables. 



Beyond pure state entanglement for atomic ensembles 



2 



The interaction between a polarized atomic ensemble and an off-resonant linearly 
polarized laser beam results - at the classical level - in a Faraday rotation. The latter 
refers to the rotation experienced by the polarization of light when it propagates through 
a magnetic medium, in our case the polarized atomic ensemble. Furthermore, this 
interaction leads to a quantum interface, i.e., an exchange of quantum fluctuations 
between the matter and the light. Seminal results exploiting such interface are the 
generation of atomic spin squeezed states produced by the interaction of coherent 
atomic sample and a squeezed light beam [HE], the realization of the atomic quantum 
memory [3] and the establishing of entangled state of two spatially separated atomic 
ensembles [I]. The above situations require a final projective measurement on the light 
(homodyne detection) that projects the atomic ensembles into their desired final state. 
This "measurement induced" mapping of fluctuations between different physical systems 
provides a powerful tool to design quantum correlations. The potential applications 
of such a genuine quantum tool have just started. Among them are the application 
in quantum metrology [ElE], quantum magnetometry [7] and quantum spectroscopy 
for detection of magnetic ordering in spinor Bose-Einstein condensates [8], or strongly 
correlated spin systems realized with ultracold atomic gases [9HT2]. 

A milestone achievement using quantum matter-light interfaces was the generation 
of entanglement between two spatially separated atomic ensembles [I] mediated by a 
single light beam propagating sequentially through both of them. In such setup, the 
entanglement between the atomic ensembles was established as soon as the light was 
measured, independently of the outcome of the measurement. However the detection 
of the entanglement thus generated, based on variance inequalities criteria [T3JH3], is 
rather challenging. As shown in jl] its verification can be made feasible with the help 
of external magnetic fields applied to the ensembles. Nonetheless, it has been shown 
in [TH] that entanglement verification without external magnetic fields is also possible if 
the light crosses through each sample at a given angle <Xj. In this setup, that we name 
"geometrical", the incident angles play a critical role since they can be used to shape 
both quantitatively and qualitatively, the entanglement between the atomic ensembles. 

The geometrical scheme of entangling different atomic systems using quantum 
interfaces can be extended beyond the bipartite pure state setting to include both 
several parties and noise. Multipartite mixed states entanglement, has been mainly 
addressed for finite dimensional Hilbert spaces and it is poorly understood in all but 
the simplest situations. For instance, recent results suggesting the potential supremacy 
of noisy quantum computation versus pure state multipartite quantum computation 
[T614T8] are highly intriguing. For continuous variables even less is known. Atomic- 
ensemble systems belong to the important class of continuous-variable states known as 
Gaussian, whose characterization and mathematical description is much simpler. For 
instance, the complete classification of three-mode Gaussian systems has been provided 
in [121120] • However, even in the simplest multipartite Gaussian scenario, the verification 
of the entanglement using variance inequalities [H] becomes a very intricate task in the 
atomic-ensemble setup. Here we address the generation and detection of multipartite 
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entanglement in Gaussian mixed states in experimentally feasible systems using the 
geometrical scheme. In particular, we provide physical bounds on atomic ensembles' 
parameters that lead to generation of bound entanglement, providing a guideline on 
their experimental accessibility. Bound entanglement refers to the entanglement which 
cannot be distilled, i.e., no pure state entanglement can be obtained from it by means 
of local operations and classical communication (LOCC). In the multipartite setting, a 
state can be bound entangled with respect to certain parties while distillable with respect 
to others. Although bound entanglement may seem useless for quantum information 
processing, it has been shown that for discrete systems it can be used for certain tasks 
(e.g. channel discrimination) or activated [2TH26] . 

A particular example of bound entanglement we are going to deal with is the so- 
called Smolin state [27] . In the discrete case, this state corresponds to a four qubit mixed 
state formed by an equal-weight mixture of Bell states that has the property that its 
entanglement can be unlocked. Moreover, it has been shown that it violates maximally 
a Bell inequality, being at the same time useless for secrete key distribution [23]. The 
Smolin state has been very recently addressed mathematically within the continuous 
variables stabilizer formalism [28J. Here we will perform its covariance matrix analysis 
that gives more insight into the nature of the state and demonstrate how such type 
of states can be experimentally realized using quantum light interfaces with atomic- 
ensembles. 

Since both the atomic ensembles and the light assume a unified continuous-variable 
description and given the inherent Gaussian character of the setting we are considering, 
we use a covariance matrix formalism [29J to describe in a compact way the generation 
and manipulation of entanglement. Such a formalism, as we shall see, enlightens the 
possibilities offered by atomic samples as a quantum toolbox to generate many different 
types of continuous variables entangled states in an experimentally accessible way. 

The paper is organized as follows. In Section |2] we review, for completeness, the 
description of matter and light as continuous variables systems, and describe the matter- 
light interaction in terms of an effective Hamiltonian and the corresponding evolution 
equations. Then we derive explicitly all the steps of the matter-light interface in terms of 
covariance matrices, symplectic transformations and Gaussian operations. In Section [3] 
we analyze the properties of the interface leading to mixed state entanglement with 
atomic ensembles. First, we consider the bipartite case in Section 13.11 setting the 
bounds on system parameters that produce mixed state entanglement. In Section I3~2l we 
move to multipartite bound entanglement in continuous variables. Here, we consider a 
general frame in which initially the atomic samples are in a thermal state with different 
amounts of noise. This allows us to set physical bounds on the generation of such type 
of entanglement in the CV scenario. In Section I3T31 we analyze the much more involved 
Smolin state [28] . We show that the geometrical setup we propose provides all the tools 
needed to generate and unlock the Smolin state, including a (quantum) random number 
generator. Finally, in Section H] we present our conclusions. 
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2. Continuous Variables systems 

2.1. The Faraday Interaction 

The setting we are considering consists of several atomic ensembles and light beams, the 
latter playing the role of information carriers between the atomic samples. At a time, 
only a single light beam interacts with the atomic ensembles. In previous proposals, 
the light beam after interacting with the atomic ensembles was measured, inducing 
entanglement between the different samples (see e.g. jHEL5l[29l[30]). Here, we will have a 
closer look at the situation in which the light beam is not measured after the interaction. 
Such action, as we will see, acts as a truly quantum Gaussian random number generator. 
We will show that through such procedure one can generate multipartite mixed states, 
in particular bound entangled states. 

As pointed out in the introduction, each atomic ensemble is described by its 
collective angular momentum J = (j x , J y , J z ). Atoms are assumed to be all polarized 
along the x direction (e.g. prepared in a particular hyperfme state) so that fluctuations 
in the J x component of the collective spin are very low and this variable can be treated 
as a classical number j x w (j x ) = h~J x = frN at j. By appropriate normalization the 
orthogonal spin components are made to fulfill the canonical commutation relation, 
jy/ y/hJ x , j z / y/hJ x = ih. Notice that they have non-zero fluctuations. To stress the 
continuous variable character of the system, we rename the above variables as "position" 
and "momentum" : 

xa = , , Pa = -j (1) 

yhJ x yhJ x 

From now on we will use only the canonical variables xa,Pa to refer to the atomic 

sample, where the subindex A stands for atomic ensemble. Later on when dealing with 

different atomic ensembles the notation XA, n ,pA,n, when we refer to the nth atomic 

sample, will be used. 

Light is taken to be out of resonance from any relevant atomic transition and 
linearly polarized along the x-direction. We use the Stokes description s = (s x ,s y ,s z ) 
for the light polarization. The components Sk (k = x, y, z) correspond to the differences 
between the number of photons (per unit time) with x and y linear polarizations, ±7r/4 
linear polarizations and the two circular polarizations, i.e., 

h h 



The Stockes operators are well suited for the microscopic description of interaction 
with atoms, however, effectively only the following macroscopic observables will be 
relevant: = J Sfc(t)dt, where T is the duration of the light pulse. The so defined 



Beyond pure state entanglement for atomic ensembles 



5 



operators obey standard angular momentum commutation rules. The assumption of 
linear polarization along direction x allows for the approximation S x ~ (S x ) = N p hh/2. 
Once more, the remaining orthogonal components S y and S z are appropriately rescaled 
in order to make them fulfill the canonical commutation rule, S y / y/HS x , S z / y/frS x 
ih. Straightforwardly, an equivalent equation to equation ([Q) arises: 

Sy A S z , . 

XL = -J=, PL = -7=, (3) 

which allows to treat the light polarization degrees of freedom on the same footing as 
the atomic variables. 

In the situation in which a light beam propagates in the Y Z plain and passes 
through a single ensemble at angle a with respect to direction z, the atom-light 
interaction can be approximated to the following QND effective Hamiltonian (see [31] 
and references therein for a detailed derivation): 

Hint( a ) = ~ k Pl{pa cos a + xa sin a). (4) 

The parameter k is the coupling constant with the dimension of the inverse of an 
action. Notice that such Hamiltonian leads to a bilinear coupling between the Stokes 
operator and the collective atomic spin operators. Evolution can be calculated through 
the Heisenberg equation for the atoms and using Maxwell-Bloch equation for light, 
neglecting retardation effects. The variables characterizing both systems (atom and 
light) transform according to the following equations ( [31] and references therein): 

x°f = x i X- nfl cos a, (5a) 
ff^S + ^sina, (56) 
x L ut = x i £-K(p i %cosa + x i A 1 sina), (5 c) 

where the subscript in/out refers to the variable before/after the interaction. The 
above equations can be generalized to the case in which a single light beam (xl,pl) 
propagates through many samples shining at the nth sample at a certain angle a n . A 
complete description taking into account inhomogeneities of the atomic samples or the 
light beams has been considered in [32] . 

Due to the strong polarization constraint, both the atomic ensembles and the 
light are initially Gaussian modes. Moreover, the Hamiltonian is linear in both atomic 
and light quadratures, and therefore quadratic in creation and annihilation operators. 
Such interaction Hamiltonians correspond to Gaussian transformations preserving the 
Gaussianity of the input state. These facts enable us to tackle the quantum atom-light 
interface within a covariance matrix formalism. 

2.2. The atom light interface in the covariance matrix formalism 

We start by reviewing the most basic concepts needed to describe Gaussian continuous- 
variable systems. For further reading, the reader is referred to [33ti35] and references 
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therein. For a general quantum system of iV pairs of canonical degrees of freedom 
("position" and "momentum"), the commutation relations fulfilled by the canonical 
coordinates R = (xi,p±, . . . , xn,Pn) can be represented in a matrix form by the 
symplectic matrix J N : [R i: Rj} = ih(J N )ij, i,j = l,..., 2N, where 

j» = ®j. j=(_°! jY («) 

Gaussian states are, by definition, fully described by the first and second moments of the 
canonical coordinates. Hence, rather than describing them by their infinite-dimensional 
density matrix g, one can use the Wigner function representation 

W(C) = ^^exp [-(C - cOV^C " d)} , (7) 

which is a function of the first moments through the displacement vector d, and of the 
second moments through the covariance matrix 7, defined as: 

di = Tr(gRi), 7^ = Tr(g{Ri - d it Rj - dj}). (8) 

The variable £ = . . . , xn,Pn) is a real phase space vector with probability 

distribution given by the Wigner function. The covariance matrix corresponding to 
a quantum state must fulfill the positivity condition 

l + iJ N > 0. (9) 

In the particular case of a physical system consisting of several atomic ensembles and 
single light beam the most general covariance matrix takes the form 

(10) 

where the submatrix j L corresponds to the light mode, 7 A to the atomic ensembles, 
that initially reads 7^ = 7^ © • • • © 7 in n and C accounts for the correlations between 
the atomic ensembles and the light. 

If a Gaussian state undergoes a unitary evolution preserving its Gaussian character, 
as it is the case here, then the corresponding transformation at the level of the covariance 
matrix is represented by a symplectic matrix S acting as 

7out = S T j in S. (11) 

Let us illustrate how to reconstruct the evolution of the covariance matrix from the 
propagation equations f )5fflj) -fl5d |) . Notice that the variables describing the system after 
interaction (out) are expressed as a linear combination of the initial ones (in). Let us 
denote this linear transformation by K 

K : (£5£, p%, xT, pTf = K(xl n , p% n , X?, p?) 7 '. (12) 
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In our case, K can be straightforwardly obtained from the evolution equations (I5aj) -(l5 
For a single atomic mode the transformation reads: 



"out 

Pa 

"out 



V PT 



J 



( 



\ 



1 

o 

— k sin a 






1 

—k cos a 




— k cos a 
ft sin a 

1 



\ 



/ 



( x 'a\ 

"in 

Pa 



(13) 



\ PL J 

Since the interaction Hamiltonian is bilinear, the matrix K can be directly applied to a 
phase space vector ( and correspondingly to the covariance matrix, however the sign of 
the coupling constant k should be changed. This is so because the phase space variables 
evolve according to the Schrodinger picture, whereas the quadratures, being operators 
transform according to the Heisenberg picture. Therefore, we define K = K\ K u_ K y 
which we apply to the phase space vector and covariance matrix as 

CutTin'Cout = Ci^Ti^Cin = (KK-HniK^-'r'Gn = C^utCin, (14) 

leading to S = (K T )~ 1 . The above formalism has been explicitly developed for a single 
sample and a single beam, but it easily generalizes to an arbitrary number of atomic 
ensembles and light beams, as well as to different geometrical settings. 

Finally, the last ingredient essential to describe the matter-light interface at the 
level of the covariance matrix is the effect of the homodyne detection of light [33J. A 
homodyne measurement on the light quadratures acts as a Gaussian map on the atomic 
covariance matrix. Assuming a zero initial displacement and covariance matrix of the 
form (FIU|) . the measurement of the quadrature xl with outcome xl leaves the atomic 
system in a state described by a covariance matrix [29 |l33"l l36] . 

</' = 1 A - CiX^X)" 1 ^, (15) 

and displacement 

d A = C(X 1 L X)- 1 (x L ,0), (16) 

where the inverse is understood as an inverse on the support whenever the matrix is 
not of full rank and X is a two-dimensional diagonal matrix with diagonal entries (1, 0). 
The measurement of the quadrature pi affects the system analogously. Let us note 
here that if the light is not measured after the interaction, the state of the atomic 
sample is characterized by the covariance matrix •y A obtained after tracing out the light 
degrees of freedom, with the displacement being a Gaussian random variable according 
to equations fl5a|) and (l5B 

d A = — k-Pl cos a, KpLsma, (17) 

where pi denotes a Gaussian random variable associated with momentum operator. 

An important step in the matter-light interface when several atomic ensembles 
are present is the analysis of the quantum correlations between the different atomic 
samples once the light beam has been measured. The symplectic formalism provides 
the whole information about the atomic covariance matrix and displacement vector after 



Beyond pure state entanglement for atomic ensembles 



8 



the interaction. This makes verification of entanglement amenable to covariance matrix 
entanglement criteria (see also [29,30j). 

An operational separability criterion, i.e., state-independent, which can be only 
applied when the full covariance matrix is available, is the positive partial transposition 
(PPT) criterion [371 [38] • F° r continuous variable systems, it corresponds to partial 
time reversal of the covariance matrix [39], i.e. a change of the sign of the momentum 
for the chosen modes. If the partially time reversed covariance matrix does not fulfill 
the positivity condition fl9]), the corresponding state is entangled. This test, however, 
checks only bipartite entanglement. For Gaussian states this criterion is necessary and 
sufficient to detect entanglement in all partitions of 1 x N modes. In the multipartite 
scenario, a state may be PPT with respect to all its bipartite divisions and still not 
be fully separable. Such states are bound entangled states. For three-mode Gaussian 
states, a operational separability criterion distinguishing a fully separable state from a 
fully PPT entangled state was given in [20] • We will use such criterion together with the 



PPT criterion in section 13.21 to demonstrate how different types of entanglement arise 
in tripartite cluster-like states at finite temperature. 

Experimentally, however, it is more convenient to check separability via variances of 
the collective observables, originally proposed for two mode states in [13] and generalized 
to many-mode states in [14] . Such criteria states that if an N mode state is separable, 
then the sum of the variances of the following operators: 

u = hi%i + . . . + h N x N 

v = gm + ■■■ + g N pN 

is bounded from below by a function of the coefficients hi, . . . , h^, gi, 
Mathematically, the inequality is expressed as 



(18) 

• 9n- 



(Am) 2 + (Av) 2 > f(h 



h 



N, 9l, 



where 



f(hi 



h>N, 9u---i 9n) 



hgi + ^2 h r9r 



rel 



+ 



h m g m + h s g t 



ser 



(19) 



(20) 



In the above formula the two modes, I and m, are distinguished and the remaining ones 
are grouped into two disjoint sets / and V . The criterion (Tl9l) holds for all bipartite 
splittings of a state defined by the sets of indices {1} U / and {m} U I'. For two mode 
states, the criterion becomes a necessary and sufficient entanglement test, however only 
after the state is transformed into its standard form by local operations [13J. This local 
transformations, however, are determined by the form of the covariance matrix. In this 
sense, the knowledge of the full covariance matrix is essential in order to determine 
whether the state is entangled. Since in experiments usually one does not have access to 
the full covariance matrix, one cannot assume that this criterion decides unambiguously 
about separability. 
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Here we analyze mixed state entanglement for atomic ensembles. We start our analysis 
with bipartite states and show that the entanglement induced by the measurement of 
light, despite its irreversible nature, can be erased by making the samples interact with 
a second light beam, in a similar fashion as it happens for the pure states [15] . We then 
address multipartite entanglement and show how bound entanglement can be created 
in a tripartite setting (cluster state) using thermal states. Finally, we analyze the effect 
of randomness introduced in the multipartite setting by the action of the light which 
interacts with all the atomic samples but is not measured. We show that through such 
procedure one may produce unlockable bound entanglement. 



3.1. Bipartite entanglement of thermal states. 

Let us start with the setup in which the atomic samples are not in the minimum- 
fluctuation coherent state (vacuum), but in a general thermal state. Under such 
assumptions, the initial state of the composite system is given by the following covariance 
matrix for atoms and light 

7 in = n 1 i 2 4 ©---©n iV i 2 4 ©i2, (21) 

where the identity I2 stands for a single mode and parameters ni, . . . , are related to 
temperature through n« = 1/ tanh[/3jU;/2] (i = 1, • • • , N), where (3 is the inverse of the 
temperature, and u is the effective frequency of the single sample. 




c) z 

Figure 1. The sketch of the setups using the geometrical approach to generate 
and manipulate bipartite entanglement. The interaction between the light beam and 
atomic samples followed by the measurement introduces squeezing in a) pa,i + Pa,2 b) 
xaa - xa,2 and c) XAA + XA.2- 

The QND Hamiltonain (j3J) for the many-mode setup reduces to 

#hrt(«) ^ ^jPL [j>A,i COSOj + X A ,i shlC^). (22) 
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where ol\ refers to the incident angle at which the light impinges the atomic sample i. 
The corresponding symplectic matrix describing this interaction is given by: 



int 



/ 




Gi 














\ Mi 


■■■ M N 





(23) 



with 



Gi 



— Ki Sin QJj 

—Ki cos a,- 




Mi 



-Kj COS Ctj 



Kj sin otj 



(24) 



For the two mode atomic states, the covariance matrix of the atomic samples after the 
interaction (setup in figure [Tb) can be straightforwardly calculated from equation fill I) 



7out 



/ ni + k 2 





K 2 








«\ 





m 













K 2 





n 2 + k 2 








K 











n 2 


n 2 n 













n 2 K 


1 + riiK 2 + ^k 2 





V « 





K 








1 / 



(25) 



Due to the atom-light interaction, both atomic modes are entangled with light, however 
the reduced state of the two ensembles is separable as one can easily check by applying 
the PPT criterion to the covariance matrix in the upper-left block. Entanglement 
between atomic samples is not produced until one measures a quadrature of light. 
Assuming the measurement outcome on xi to be xl,i> the covariance matrix describing 
the final state of the samples is given by [see (fl5j) and ffT6]) ] 



Tfin 



/ Hi + K 2 



2 



V 



K 







nin2K 2 +ni 
(rai+ra2)K 2 +l 



(ni+ri2)« 2 +l 



K 



n 2 + k 2 








(n\+n2)K 2 +l 



nin 2 K 2 +ri2 
(ni+ri2)K 2 +l 



(26) 



and the displacement of the final state is 



d 



fin 



0. 



,0, 



(27) 

(m + r^ + l'"' { ni +n 2 )K 2 + 1* * 
For what follows it is important to notice that the covariance matrix is independent of 
the measurement outcome, but the latter is clearly present in the displacement vector 
of the atomic modes. To check the entanglement between the atomic samples after 
the light has been measured we use the separability criterion based on the variances of 
the two commuting operators [13J, which states that for any separable state the total 
variances fulfill 



-i 2 



A(|A|zU,i + jPa,2) 



+ 



1 2 



A(\X\x Ay 



A 



X A,2, 



> 2h. 



(28) 



This is a sufficient but not necessary condition for separability. 
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We restrict our analysis to a single inequality involving the collective observables 
with A = 1 in (|28p since it is the one applicable experimentally [3]. The way to measure 
such combination of variances has been described in detail in [T5] . 

Extracting the variances from the elements of the final covariance matrix §3o$) we 
obtain 

1 2 1 n l + n 2 

IA(PA,1 + PA,2>\ = o (7fin,22 + 7fin,44 + 27 firli 24) 



ft 2 2 (ni + n 2 ) K -\- 2 

1 2 1 1 

- [A(x A ,i - x A , 2 )] = TjWfiMi + 7fin,33 - 27 fin ,i 3 ) = 2 ( n i + n 2j • (29) 

The substitution of the above expressions in equation fl28l) leads to the violation of the 
separability criterion for the values of ni,n 2 , and k fulfilling 

9 2(ni + n 2 -2) , , 

(4 - m - n 2 J (ni + n 2 J 

, 2(ni + n 2 -2) , , 

k< Ta — h L V ™i+™2>4 31 

(4 - m - n 2 J (ni + n 2 ) 

One immediately notes that for n± + n 2 > 4, the inequality can never be violated, since 

the right-hand side of the inequality becomes negative. In Figure [2] we compare the 

complement of the set defined by f l3"Uj) (states that are not detected) with the set of 

separable (PPT) states. This shows that the variance inequality does not detect all 

entangled states. This is due to the fact that only one combination of variables, i.e., 

Pa,i + Pa,2 is squeezed. In this case the thermal fluctuations still present in xa,\ — xa,2 

does not allow the sum in equation f[2"51) to violate the bound. 

In the case of initially pure states an arbitrary coupling produces entanglement [T5] . 
Here we see that for initially thermal states this is not the case and the generation of 
entanglement requires a stronger coupling with light. 

In order to increase the amount of entanglement and improve its detectability 
through the variance inequality criterion, a sequence of steps displayed in Figure 
and [lb is required. These introduce squeezing in two commuting combinations of 
quadratures. Calculation similar to those from the previous paragraph lead to 

5 |A( ^' + ^ = I |A(i '« - i « )|2 = 2( Bl T^ + 2 - (32) 
and the violation of the variance's inequality ( |28|) for A = 1 is now obtained for a much 
lower coupling: 

2 wi + n 2 - 2 

k > . 33 

2n x + 2n 2 y 1 

The set of states detected by the spin variance inequality is compared in figure [2b to 

the set of separable (equivalently PPT) states. Again the spin variance inequality does 

not detect all entangled states, however, now is much more efficient than in the previous 

case since the fluctuations in both combinations of variables were suppressed. 

Interesting enough such measurement induced entanglement between two atomic 

samples can be deleted by exploiting the squeezing and anti-squeezing effects produced 

by the laser beams for pure state entanglement [15]. Here we demonstrate that this 
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Figure 2. Comparison of the sets of parameters n,n\,n2 for which the state of two 
ensembles remains separable (brown inner region), with the ranges of parameters for 
which it is not detected by the spin variance inequality (|28[) with A = 1 (gray outer 
region). In figure a) we consider the state produced in the setup in figure [T^,, whereas 
in figures b) the one produced in two steps schematically depicted in figures [TJl and 

Eb. 



is also true when the initial states are thermal, however, the procedure has to be 
slightly modified. Let us consider the state represented by the covariance matrix f l2"B"j) . 
Interaction with a second light beam impinging on each atomic sample at ot\ = a 2 = vr/2 
(as depicted in Figure [lb) followed by the light measurement will erase the entanglement 
produced by the interaction with the first light beam if the properties of light are 
appropriately adjusted. To this aim, it is sufficient to impose a second light beam 
characterized by the covariance matrix 

:ax x ) 2 





1L 




with 



(Axx) 2 
(A Pl ) 2 



K 2 n\ + K 2 n 2 + n\n 2 
n x n 2 



K 2 rii + K 2 n 2 + 1 

and interaction coupling k. The covariance matrix of the final state is given by: 



f K 2 n 2 +?ii(K 2 +n 2 ) 



7e 



V 



2K 2 +n 2 









ni(2K 2 n 2 +l) 
K 2 (ni+ri2)+l 
















\ 



K 2 +n 2 ) 



and the corresponding displacement is: 



cL 



2ft 2 +ni 





X L oK 





n 2 (2K 2 ni+l) 
K 2 (ni+n 2 )+l / 



(34) 



(35) 



(36) 



2k, 2 + n 2 ' k, 2 (ni + n 2 ) + 1 ' 2k, 2 + rii ' k 2 (ni + n 2 ) + 1 



(37) 
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where xl,2 is the measurement outcome. Note that the state fl36l) is separable. It will 
be identical to the initial one, however, only if n 2 = n x . 

3.2. Multipartite bound entanglement 

We move now to the truly multipartite entanglement and focus on the simplest case 
of the 3 atomic samples. We aim at analyzing bound entanglement, which exist 
only in the mixed state case. The classification of tripartite entanglement in CV was 
given by Giedke and coworkers in [20], where five classes of states were distinguished. 
Following their classification, inseparable states with respect to every bipartite splitting 
are denoted as class 1. States which are biseparable with respect to one (and only one) 
bipartition belong to class 2. States that are biseparable with respect to two or three 
bipartitions, but still entangled, belong to class 3 and class 4, respectively. These two 
classes of states are bound entangled. Finally, class 5 are the separable states. 

Continuous-variable cluster-like states can be a universal resource for optical 
quantum computation [4TH-42J . They are defined by analogy with the discrete cluster 
states generated via Ising interactions between qubits [32] • One can associate the modes 
of the TV-mode (CV) system with the vertices of a graph G. And the cluster corresponds 
to a connected graph [33]. For continuous variables, cluster states are defined only 
asymptotically as those with infinite squeezing in the variables 



for all the modes belonging to the graph, where N a denotes the set of neighbors of vertex 
a. Cluster-like states are defined when the squeezing is finite. 

It is possible to create such states with atomic ensembles in the proper geometrical 
setup [15] taking the samples in the initial vacuum state. The protocol consists of 
interaction with light passing through the samples at specified angles (see Table [Q plus 
a homodyne detection of light. Here we check the properties of the states obtained 
through the same procedure, however for initially thermal samples, i.e., 7^ = ni e © i 2 . 
We demonstrate that in such setup bound entangled states are produced for certain 
choices of temperature and coupling. The presence of undistillable entanglement in 
thermal finite-dimensional systems was considered in [44] . 

We analyze the two cluster states, the linear one and the triangular one, focusing 
on how the temperature destroys the genuine tripartite entanglement. To discriminate 
class 1, 2 and 3, the PPT criterion is enough. To discriminate between class 4 and 5 
we will use the operational necessary and sufficient criterion for full separability given 



The results we obtain are summarize in Figure [3] where we depict as a function 
of the temperature T through n(T) = 1/ tanh (u/2T) and coupling k, the different 
entanglement types created between the atomic samples. 

For the linear cluster state (figure [3^), for a fixed value of the parameter k, the 
state is NPT with respect to all the three cuts in the interval < T < T a , meaning that 




(38) 



beN a 



in [20]. 
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Table 1. Parameters for tripartite cluster states and the four-partite Smolin state. 
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Figure 3. Range of parameters for (a) the linear and (b) the triangular cluster state 
for which it belongs to different entanglement classes (see explanation it the text). 



it belongs to class 1. For T a < T < T\, the state is PPT with respect to the bipartition 
2 — (1,3) and NPT with respect to others. Hence, in it belongs to bound entangled 
class 3. For T > Tb the state becomes PPT with respect to all the cuts. Using the 
separability criterion found in [20J, we find that within the range < T < T c the state 
is class-4, while for T > T c it becomes fully separable. 

Figure [3b corresponds to the different entanglement classes for a triangular cluster 
state. In this case, because of the permutational invariance, class 3 is never recovered. 
Nevertheless a class-4 region still appears between the genuine tripartite and the fully 
separable states. 

3.3. Multipartite unlockable bound entanglement 

Here we propose to use the Gaussian randomness that is inherently present in the system 
due to uncertainty of the outcome of the light measurement. We will show that it is 
essential to generate the bound entangled Smolin state. 

The Smolin state introduced in [27| using qubits is an interesting example of 
an undistillable multipartite entangled state that can be unlocked if several parties 
perform a collective measurement and send the outcome to the remaining parties. The 
Smolin state if form by an equal-weight mixture of products of the four Bell states 
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^(i = l,...,4) 

4 

femolin = ^ I^XV'iUfi ® IV'iXV'tlciJ- (39) 
t=l 

One sees that parties AB and CD share the same Bell state but none of the pairs 
know which one it is. Such state is clearly separable with respect to partition AB\CD. 
Moreover, since it is permutationally invariant, the above statement also holds for every 
two pairs, not only for AB and CD. Hence, according to the definition, the state is 
undistillable with respect to arbitrary two parties, since there always exists a separable 
bipartition dividing them. 

The fact that arbitrary two pairs share the same unknown Bell state makes the 
entanglement present in the Smolin state unlockable. Imagine the situation in which 
two of the parties, say C and D, meet in one laboratory and perform collectively a 
measurement in the Bell basis, identifying thus the state they possess. This measurement 
projects the state shared by AB on the one corresponding to the outcome of the 
parties CD. Hence, if A and B, still being in separate laboratories, obtain (by classical 
communication) the information about the outcome of CD, they can use their maximally 
entangled state. We emphasize that the otherwise unknown Bell state is useless. 

The Smolin state possesses many other interesting features. Apart from being 
unlockable, its entanglement can be activated [15] through a cooperation of five parties 
sharing two copies of the state. Despite being bound entangled, this state can be used 
in the remote quantum information concentration protocol [16] and maximally violates 
a Bell inequality being at the same time useless for the secret key distillation [23]. 

Recently, Zhang has generalized the discrete Smolin state to the Gaussian CV 
scenario using the mathematical formalism of stabilizers [28]. It turns out that the 
state proposed by Zhang can indeed be understood as some mixture of products of the 
continuous variables EPR-like states in this way. Let us consider two EPR(-like) states 
with known reference displacement and randomly displace both of them in the opposite 
unknown directions in the phase space, as schematically depicted in figure HI The only 
constraint imposed on the random displacement is that it has a Gaussian distribution. 
In this way we obtain a Gaussian mixture of two displaced EPR pairs. 

The properties of the CV Smolin state are slightly different from those of the 
qubit state. Firstly and most importantly, we loose the permutational invariance [28J. 
Secondly, since the maximally entangled EPR-state corresponds to infinite squeezing in 
the CV and is not physical, we will always have to deal with "EPR-like" states and 
therefore it is necessary to analyze the effect of the finite squeezing on the properties of 
the final Smolin-like state. Both of the above constraints force the Smolin-like state to 
be undistillable only with respect to chosen parties and within a particular range of the 
squeezing parameter. 

Let us demonstrate how the mathematical setup generating the Smolin-like state 
proposed in [SB] greatly simplifies in the atomic-ensemble realization. We characterize 
each step of the procedure using the covariance matrix and the displacement vector. 
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Figure 4. Schematically depicted Smolin state in the phase space. The red dots 
correspond to the EPR-like states, the green ellipses represent the random Gaussian 
displacement that depends on the properties of light passing through the samples. 

Finally, using the symplectic formalism, we demonstrate the unlocking protocol. 

From now on we will distinguish the parties (modes) with Latin numbers. As 
demonstrated by Zhang, the Smolin-like state can be generated in two steps: 

(i) generate two EPR-like pairs with reduced fluctuations (squeezing) in x\ + x 2 and 
Pi ~ P2 (£3 + £4 and P3 — p4 for modes 3 and 4) and known displacement 

(ii) displace the EPR pairs by random Gaussian variables Ax and A2 such that the 
quadratures transform in the following way: 

xi + x 2 + £3 + £4 -> (x\ - Ai) + (x 2 - Ai) + (x 3 + Ai) + (£4 + Ai), (40) 

and 

Pi - P2 + P3 - Pi -> (pi + A 2 ) - {P2 - A 2 ) + (pa - A 2 ) - (P4 + A 2 ). (41) 

Note that indeed the squeezed combinations in modes 1 — 2 and 3 — 4 are displaced in 
the opposite directions. Implementation of each of these steps is experimentally feasible 
in the geometrical setup of atomic ensembles. Step (i) concerns the generation of EPR- 
like states with known displacement which is easily achieved (see also [15]). Step (ii) 
concerns the generation of a random Gaussian variables in the displacements. In our 
geometrical setup this can be achieved by shining two light beams at some given angles 
as summarized in table [TJ A straightforward application shows that the angles reported 
in the above table modify the the quadratures of atomic spin ensembles in the following 
way: 

X1+X2+X3+X4 -> (x 1 -Kp LA ) + (x 2 -Kp LA ) + (x 3 + Kp LA ) + (x i + Kp LA )(4:2) 

and 

P1-P2+P3-P4 -»■ (Pl+KPL,2)-{fa-KPL,2) + (fa-KPL,2)-{fa + KPL,2)-{43) 

Comparing now equations fj4"0"j) and (14T1) to equations (jl2l) and (14"3"|) respectively, one 
observes that our setup reproduces precisely the required displacements. Now, if the 
output light is not measured, npL,i and npi,2 remain unknown random Gaussian variables 
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with expectation value and variance that can be easily adjusted by choosing the input 
light with specific properties (vacuum squeezed, thermal etc.). From now on we will 
denote this random variables by Pl,i,Pl,2- 

We switch to the covariance matrix formalism that allows us to get a better insight 
into the process and the properties of the final state. We begin directly with two EPR 
pairs with zero reference displacement and characterized by the squeezing parameter r, 
and a general light mode characterized by the variances: (A£ £1 ) 2 , (Ap L ^) 2 and zero 
displacement. The initial state of the setup is described by: 

\2 



7i m) = TepRi © 7epr 2 © 



(A*L,ir 









(APL,1 



(44) 



The interaction between the atomic ensembles and the light leads to the following CM 
for atoms 
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(45) 



with 
ai = 



k 2 (Ap Ljl ) 2 + ch2r, b = ch2r c\ = k 2 (Ap Ljl ) 2 — sh2r, 
d = k 2 (Ap Lj i) 2 - ch2r e x = -k 2 {Ap LA ) 2 . 



Interaction with a second mode of light characterized by variances (Ax Lt2 ) 2 , (Ap Lt2 ) 2 , 
leads to the CM: 
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(46) 



ci = (Ap Lt i) 2 — sh2r 
ei = -k 2 (Ap LA ) 2 , e 2 = /t 2 (Ap L;2 ) 

and the displacement of the final state is given by: 

,(out) 



c 2 = k 2 (Ap Lj2 ) - sh2r, 

2/ a- \2 



4 UULJ = 



K ^>2, -kPl,i, -KPL,2, KPL,1, -KPL,2, KPL,1, KPL,2) ■ (47) 
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We verify separability of the obtained state using the partial time reversal criterion and 
see that the state: 

(i) it is always PPT with respect to partition 12 1 34, 

(ii) it is always NPT with respect to partition 13 1 24, 

(iii) it is PPT with respect to partition 14|23 only if 



k 2 (Ap 



L,lj 



e- 2r ) or k 2 (Ap Li2 ) 



2r 



-It 



); 



(48) 



(iv) all the partitions one vs. three modes are always NPT, in agreement with the 
results obtained in [28J. 

An important thing to notice is that the state is unlockable bound entangled only if the 
partition 14 1 23 is PPT, since only in this case the entanglement between modes 1 and 
2 is bound and the unlocking procedure makes sense. 

In order to show the unlockability of the obtained state, we write the 
CM and displacement in the basis in which the EPR states are diagonal, i.e., 
(l/\/2) (xi + x 2 , ■ ■ ■ ,p3 — Pa). For simplicity we also assume that (Ax L ^) 2 = (Axl j2 ) 2 = 
{Ax L ) 2 and (Aj3 Lj i) 2 = (Ap Lj2 ) 2 = {Ap L ) 2 then 
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(49) 



+ fj 



(50) 



Now inspection of 7 S shows clearly that we have two correlated EPR pairs displaced by 
random vectors pointing in the opposite directions in the phase space of the squeezed 
variables, exactly as depicted in figure HI Since we do not know where in the phase 
space are the EPR pairs, we cannot use them for quantum tasks. Therefore, by analogy 
to the discrete case [221EZ], unlocking of entanglement is understood as learning which 
of the displaced EPR pairs is shared by pair 1 — 2 by measuring the displacement of 
the other EPR pair shared by 3 — 4. 

For completeness we also demonstrate all the steps of the unlocking protocol with 
the atom-light interface using the symplectic formalism. The measurement of X3 + £4 
and ps —pi can be done with two probe light beams with reduced fluctuations in xp, i.e., 
(Axp) 2 < h/2 (QND measurement) jl?]. The probe beams pass through the samples 



3 — 4 at the angles a 
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7T, respectively, and are 



measured afterwards. If the measurement outcomes are, respectively, x + , p-, the final 
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state of modes 1 — 2 is characterized by the following CM and displacement (still in the 
EPR basis): 

f ,f(2^ + ^(AM 2 ) 2r n \ 

e 2r (2K 2 f+(Ax P f )+2K 2 ^ 
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(52) 



(Ax P ) z + 2K 2 (e~ 2r + f) 

If the variables x^ + £4 and j3 3 — p 4 are strongly enough squeezed, the obtained outcome 
approximates well the displacement of X\ + x 2 and p\ — j3 2 . Therefore, we may write 
2^Pl,i = — and 2Kpi^ = —p~- In this way the resulting displacement is a function 
of the known parameters 
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(53) 



In comparison to the initial EPR state, the covariance matrix (15T|) has an additional 
positive term in diagonal elements, which reads 

/ (2k 2 + e 2r (Axp) 2 ) 



6 



(54) 



e 2r (2k 2 f + (Axp) 2 ) +2k 2 ' 

It increases the variance of the squeezed variables and therefore makes the state mixed 
and less entangled. In order to quantify the entanglement in the final unlocked state, 
we assume that the coupling constants k used in the preparation of the EPR states and 
the ones used for generation and unlocking of the Smolin state are equal and we express 
them in terms of the squeezing parameter r 



K 



[l + e 2r )/2. 



(55) 



Further we assume that the squeezing of the probe light beam is the same as the 
squeezing of initial EPR pairs, i.e., (Axp) 2 = e~ 2r . In this way, the properties of 
the final state depend only on two parameters: r and (Aj3l) 2 . In figure [5] we plot the 
negativity of the unlocked state, as a function of these parameters. For the admissible 
values of r, the negativity of the final state is roughly a half of the value for the initial 
EPR pair. 



4. Conclusions 



Our results can be summarize us follows. We have shown that, unlike in the pure state 
case, if the atomic ensembles are initially in a thermal state the entangling procedure 
only succeeds for certain choices of the coupling parameter. In this way we have provided 
bounds on the experimental parameters leading to mixed state entanglement. Further 
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Figure 5. Negativity of the unlocked state given by equation ((52)) as a function of 
the squeezing parameter r and variance of the random displacement (Ap^) . All the 
other parameters are expressed in terms of r (see explanation in the text). 

we have shown that tripartite continuous cluster-like states are robust against initial 
thermal noise present in the atomic ensembles, meaning that there exists a finite range 
of temperatures for which the final state remains genuine tripartite entangled. We have 
also demonstrate that with increasing temperature, the states becomes bound entangled 
and finally fully separable. 

Finally, we have also shown that a mixed state of atomic ensembles can be produced 
from an initial pure state if the light beam is not measured after the interaction. In this 
case, the interaction with the light beam mixes all the possible outcomes with a Gaussian 
weight. Using this procedure, starting from two pure bipartite entangled states, it is 
possible to produce the so-called Smolin state which is an example of an unlockable 
bound entangled state. It should be emphasized that, since this procedure does not 
involve the final projective measurement, it cannot produce extra entanglement. Its 
only effect is introducing randomness in the system, making the entanglement bound. 

Summarizing, our analysis shows that atomic ensembles offer a versatile toolbox for 
generation and manipulation of multipartite entanglement both in the pure and mixed 
state cases. By exploiting further the geometrical setting introduced in [15] we have 
been able to assess the robustness of the entangling schemes to set guidelines for future 
experiments, both in the bipartite as well as in the multipartite case. Moreover, we 
have demonstrated that the very same setting provides experimentally feasible schemes 
for generation of bound entanglement, broadening the applicability of atomic ensembles 
for quantum information studies. In particular, our results provide, to the best of 
our knowledge, the first experimentally feasible realization of the Smolin state with 
continuous variables (for qubit realization see |48j). 
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